Prom stress-induced fluidization processes 
to Herschel-Bulkley behaviour in simple yield stress fluids 
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Stress-induced fluidization of a simple yield stress fluid, namely a carbopol microgel, is addressed 
through extensive rheological measurements coupled to simultaneous temporally and spatially re- 
solved velocimetry. These combined measurements allow us to rule out any bulk fracture-like sce- 
nario during the fluidization process such as that suggested in [Caton et at, Rheol Acta, 2008, 
47, 601-607]. On the contrary, we observe that the transient regime from solidlike to liquidlike 
behaviour under a constant shear stress a successively involves creep deformation, total wall slip, 
and shear banding before a homogeneous steady state is reached. Interestingly, the total duration 
Tf of this fluidization process scales as tj oc 1/(<t — er c )' 3 , where a c stands for the yield stress of 
the microgel, and ft is an exponent which only depends on the microgel properties and not on 
the gap width or on the boundary conditions. Together with recent experiments under imposed 
shear rate [Divoux et al., Phys. Rev. Lett., 2010, 104, 208301], this scaling law suggests a route 
to rationalize the phenomenological Herschel-Bulkley (HB) power-law classically used to describe 
the steady-state rheology of simple yield stress fluids. In particular, we show that the steady-state 
HB exponent appears as the ratio of the two fluidization exponents extracted separately from the 
transient fluidization processes respectively under controlled shear rate and under controlled shear 
stress. 

PACS numbers: 83.60. La, 83.50.Ax, 83.50.Rp 



I. INTRODUCTION 



Yield stress fluids (YSF) are widely involved in man- 
ufactured products such as creams, gels, or shampoos. 
These materials are characterized by a transition from 
solidlike to liquidlike above the yield stress a Cl which is 
of primary importance at both the manufacturing stage 
and the end- user level [l| . Recently it was recognized that 
simple YSF, which mainly consist in emulsions, foams, 
and carbopol microgels, can be clearly distinguished from 
thixotropic YSF 2] : in steady state the former ones can 
flow homogeneously at vanishingly small shear rates un- 
der controlled stress while the latter exhibit a finite 
critical shear rate [5|, |6| . Still, in spite of its importance 
for applications, the transient fluidization process of sim- 
ple YSF has remained largely unexplored and previous 
works have focused either on global rheometry under an 
applied stress [l], 0> Hi or on time- resolved local velocime- 
try under controlled shear rate [9l4ll|. Thus, detailed 
local information concerning the fluidization of a simple 
YSF under applied shear stress are still lacking, which 
prevents to make clear connections with observations un- 
der imposed shear rate and with steady-state rheology. 

In this article we report a temporally and spatially 
resolved study of the stress-induced fluidization of car- 
bopol microgels through ultrasonic echography. Our aim 
is to address the following basic questions: (i) What is 
the fluidization scenario of such a simple YSF under im- 



posed shear stress? (ii) How does it compare to imposed 
shear rate experiments? (iii) Can one make a connec- 
tion between these transient fluidization processes and 
the steady-state rheology, which is well described by the 
Herschel-Bulkley (HB) law [|, HM1? Here, we show 
using an ultrasonic velocimetry technique that carbopol 
microgels submitted to a constant shear stress a un- 
der rough boundary conditions successively exhibit creep 
deformation, total wall slip, and shear banding before 
reaching a homogeneous steady state. A close inspection 
of the backscattered ultrasonic signals allows us to rule 
out a scenario involving bulk fracture of the material. 
The duration of the fluidization process decreases as a 
power-law with the reduced shear stress a — a c . This 
power law only depends on the sample preparation pro- 
tocol and not on boundary conditions or on the cell gap. 
Together with recent experiments under imposed shear 
rate [Tl| , this provides for the first time a direct link be- 
tween the yielding dynamics of a simple YSF and the HB 
law which accounts for its steady-state rheology. 
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FIG. 1: Top: Shear-rate 7 vs time t for a shear stress a = 37 Pa applied at t = (batch 1, rough BC). Inset: zoom on the 
shear-banding regime. The dashed curve shows d"y(t)/dt whose maximum corresponds to r/. Bottom: velocity profiles v(r), 
where r is the distance to the rotor, at different times [symbol, time (s)]: (Oi 500), (A, 3460), (V, 5890), (•, 6450), (o, 6620), 
(♦, 17349). On each profile, the upper value of the velocity scale is set to the current rotor velocity no so that the slip velocity 
v s at the rotor can be read directly as v s — vq — v(r — 0). 



II. EXPERIMENTAL 
A. Sample preparation 

1. Working system: carbopol microgels. 

Our working system is a microgel made of car- 
bopol ETD 2050 which comprises homo- and copoly- 
mers of acrylic acid highly cross-linked with a polyalkenyl 
polyether Q101. The 

microgel is traditionally prepared 
in two steps: (i) the polymer is dispersed in water leading 
to a suspension of carbopol aggregates and (ii) a neutral- 
izing agent is added (in our case sodium hydroxyde) lead- 
ing to polymer swelling and to microgel formation. The 
microstructure of such microgels consists in an assembly 
of soft jammed swollen polymer particles, with typical 
size ran ging ro ughly from a few microns to roughl y 2 
microns |16l - [l8j depending on the type of carbopol [15j . 
its concentration [13| , the final value of the pH [19| , the 
type of neutralizing agent [ljj and, last but not least, 
the stirring speed during neutralization fl5| . Carbopol 
microgels exhibit good temperat ure stability [2(|. They 
are also known in the literature to be non-aging, non- 



thixotropic simple YSF 0, [ll [HI, HI and their steady- 
state flow curve nicely follows the HB law: 

a = a c + fjj n , (1) 

where 7 is the shear rate and n = 0.3-0.6 dependin g on 
the type of carbopol and its concentration H [ll|, [ill.Tl8j . 



2. Sample preparation protocol. 

For our study we prepare two kinds of samples: tra- 
ditional samples on the one hand, and samples that are 
seeded with micronsized glass spheres on the other hand, 
in order to use ultrasonic speckle velocimetry (USV) [23[ 
simultaneously to standard rheological measurements. 
The detailed protocol for seeded samples is as follows: 
we first add 0.5% wt. of hollow glass spheres (Potters, 
Sphcriccl, mean diameter 6 /im, density 1.1) in ultrapure 
water; pH increases roughly from 7 to 8. As carbopol is 
hydrosoluble only for pH< 7, we add one or two drops 
of concentrated sulfuric acid (H2SO4, 96 %) to make the 
pH drop to roughly 5. The glass sphere suspension is 
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heated at 50° C and the carbopol powder is carefully dis- 
persed under magnetic stirring at 300 rpm for 40 min, at 
a weight fraction C, with C ranging from 0.5 to 3% wt. 
The mixture is then left at rest at room temperature for 
another 30 min, after which pH ~ 3. Finally we neu- 
tralize the solution with sodium hydroxide (NaOH, con- 
centration 10 mol.L -1 ) until pH = 7.0 ± 0.5 while stir- 
ring manually. This leads to a carbopol microgel which 
is finally centrifuged for 10 min at 2500 rpm to get rid 
of trapped bubbles. As for traditional samples without 
glass beads, the protocol starts directly by adding the 
carbopol powder to a heated volume of ultrapure water 
and continues as explained above. 



3. Influence of the preparation protocol on the batch 
properties. 

We emphasize that the final microgel macroscopic 
properties are quite sensitive to the preparation proto- 
col. In particular, it is well known in the literature that 
during the neutralization step, (i) the way the base is 
added (drop by drop or all at once) as well as (ii) the 
exact final value of the pH and (iii) the stirring speed 
during the neutralization process influence the values of 
the parameters of the HB model [IE |24j|. This derives 
from the fact that these three parameters control the 
particle size distribution of the microg el HE G!. 

In this 

paper, we take good care to neutralize our samples in a 
reproducible fashion. Nonetheless, from batch to batch, 
the final pH value of the microgel varies in the range 
6.5 < pH < 7.5. Therefore, when comparing different 
geometries, gaps, or boundary conditions, we pay spe- 
cial attention to use results obtained on a single batch, 
so that the preparation protocol does not introduce any 
bias. We will mainly report data obtained on two dif- 
ferent batches of carbopol weight fraction C = 1 % wt. 
and seeded with glass spheres, noted batch 1 and batch 
2, prepared separately but following the same protocol. 
We will also discuss the influence of the carbopol concen- 
tration C on four different traditional unseeded batches 
prepared separately: C — 0.5, 0.7, 1, and 3 % wt. 



4- Influence of the seeding glass spheres. 

Linear viscoelasticity measurements show that the ad- 
dition of hollow glass spheres generally slightly stiffens 
the system (by at most 10%)|25[. However, we shall check 
throughout the whole manuscript that traditional sam- 
ples and seeded samples exhibit the exact same rheologi- 
cal trends, which demonstrates that the acoustic contrast 
agents play no significant role in the fluidization scenario 
under imposed shear stress. 



B. Experimental setup and protocol 

1. Rheological setup. 

Rheological measurements are performed with a stress- 
controlled rheometer (Anton Paar MCR301). Two differ- 
ent small-gap Couette cells were used to test the influence 
of the boundary conditions (BC) on the fluidization pro- 
cess: a rough cell (surface roughness 5 ~ 60 /Ltm obtained 
by gluing sandpaper to both walls, rotating inner cylin- 
der radius i?i n t = 23.9 mm, gap width e = 1.1 mm, 
and height h = 28 mm) and a smooth Plexiglas cell 
[S ~ 15 nm (AFM measurements), R m t — 24 mm, 
e = 1 mm, and h = 28 mm]. Also, to test the influence 
of both the geometry and the gap, experiments were per- 
formed with a plate-plate geometry (radius 21 mm, gap 
width e = 1 and e = 3 mm) with two different bound- 
ary conditions: rough (glued sandpaper, 8 ~ 46 /im) and 
smooth [glass, S ~ 6 nm (AFM measurements)]. Finally, 
note that for both geometries, we use a solvent trap in- 
cluding a cover and a small water tank so as to efficiently 
prevent evaporation. 

2. Local velocity measurements. 

Velocity profiles are measured at about 15 mm from 
the cell bottom through ultrasonic speckle velocimetry 
(USV) as described in details by Manneville et al. [23| 
In brief, USV relies on the analysis of ultrasonic speckle 
signals that result from the interferences of the backscat- 
tered echoes of successive incident pulses of central fre- 
quency 36MHz generated by a high-frequency piezo- 
polymer transducer (Panametrics PI50-2) connected to 
a broadband pulser-receiver (Panametrics 5900PR with 
200 MHz bandwidth). The speckle signals are sent to a 
high-speed digitizer (Acqiris DP235 with 500 MHz sam- 
pling frequency) and stored on a PC for postprocessing. 
A cross-correlation algorithm yields the local displace- 
ment from one pulse to another as a function of the radial 
position across the gap with a spatial resolution of 40 /im. 
After a calibration step using a Newtonian fluid, tangen- 
tial velocity profiles are then obtained by averaging over 
10 to 1000 successive cross-correlations depending on the 
desired temporal resolution. During the creep experi- 
ments considered here, the shear rate 7 strongly varies 
in time so that the technique described in [23| for a fixed 
shear rate has to be slightly modified: we use the analog 
output of the MCR 301 rheometer to monitor the rotor 
velocity in real-time and to constantly adapt the repeti- 
tion frequency of the ultrasonic pulses to the current -y(t) 
in order to follow the dynamics. 

3. Experimental protocol. 

Before starting an experiment on a fresh sample, a 
strong preshear is applied for 1 min at +1000 s _1 and 
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FIG. 2: (a) Shear rate 7 vs time for a = 35, 36, 37, 38, 40, 
43, 45, 50, 55, 60, 65, and 70 Pa (from bottom to top) for 
batch 1 and rough BC. (b) Fluidization time 77 extracted 
from the second inflection point of j(t) (o) and from USV (•) 
vs the reduced shear stress a — a c with a c = 28.3 Pa. The 
gray line corresponds to the best fit r/ = B/(a — a c ) with 
B = (3.77 ± 0.07) ■ 10 7 and /3 = 4.0 ± 0.1. The inset shows 
the least square determination of a c (see text). 

for 1 min at -1000 s -1 to erase the loading history. The 
viscoelastic moduli are then monitored for 2 min. We 
found that both the elastic and the viscous moduli do 
not vary significantly after 2 min. Finally, the sample 
is left at rest for 1 min to ensure that a reproducible 
initial state is reached. Let us note that we used the ex- 
act same protocol for experiments performed under con- 
trolled shear rate [ll[ so as to compare results obtained 
for both applied shear stress and applied shear rate. 

III. RESULTS 

A. General scenario under rough boundary 
conditions 

Velocity profiles are measured simultaneously to the 
global shear rate j(t) during long creep experiments un- 
der an imposed shear stress a. As shown in Fig. Q] for 
(7 = 37 Pa and rough BC, the shear rate follows an S- 
shaped curve with an initial decay characterized by a 
power law j(t) oc i -2 / 3 , which is strongly reminicent of 



the Andrade creep observed in plastically deforming crys- 
tals [H, H3| and some viscoelestic materials [28] . This 
slow decay is followed by a strong increase of i(t) before 
a steady-state is reached after about 10 4 s. Local ve- 
locimetry allows us to further distinguish four regimes in 
creep experiments. The microgel first experiences creep 
deformation for t < 1500 s with extremely small veloci- 
ties (< 1 /im.s -1 ) [Fig. Ufa)]. At t ~ 1500 s, the microgel 
fails at the inner wall and, in spite of rough BC, it en- 
ters a regime of total wall slip which lasts until t ~ 4700 s 
[Fig. Hp))]. For 4700 ^ t ^ 6600 s, a shear band nucleates 
at the rotor and grows across the gap [Fig. QJc)-(d)]. Full 
fluidization of the sample is achieved at t — Tf ~ 6600 s, 
after which the velocity profiles remain linear while the 
slip velocity vanishes [Fig.[Tfe)-(f)]. Interestingly, we ob- 
serve that complete fluidization coincides almost exactly 
with the second inflection point of the curve 7(2) [inset 
of Fig. [T] so that the fluidization time r/ may also be 
defined as the time at which d7(t)/dt passes through a 
maximum. As seen in Fig. [5Jb) and Fig. (3{b), these two 
independent estimations of Tf lead to negligible discrep- 
ancies. 



B. Influence of the boundary conditions 

As reported in Fig. [3 and Fig. [31 creep experiments 
were repeated for a ranging from 35 to 70 Pa (32 to 60 
Pa resp.) for rough (smooth resp.) BC . The above flu- 
idization scenario is very robust and also holds for smooth 
BC. In this last case however, no creep deformation is ob- 
served and the slip regime starts almost instantly. Corre- 
spondingly j(t) does not show any power-law behaviour 
at short times [Fig. |3Ja)]. Nonetheless the rest of the 
process is identical to the one described for rough BC. 
Moreover, as seen in Fig. [2jb) and Fig. [3fb), the du- 
ration of the transient regime sharply decreases as the 
stress is increased. We found that the best way to fit the 
tj vs a data is a power law of the shear stress reduced 
by some well-defined critical stress o~ c : 

T f = B/{o-a c f. (2) 

The value of a c chosen here is the one that minimizes 
the x 2 °f linear fits of lnr/ vs ln(er — er c ) [see insets of 
Fig. [2Jb) and Fig.[3Jb)]. Since the sample will never flu- 
idize for a < er c , this critical stress can be interpreted as 
the yield stress of the material. At this stage, it is im- 
portant to note that our least-square procedure for mea- 
suring o~ c from the fluidization time completely differs 
from any standard way of estimating the (static) yield 
stress under controlled shear stress or from short-time 
stress overshoots under controlled shear rate [25[. We 
shall show below that this a c indeed nicely coincides with 
the (dynamic) yield stress inferred from steady-state flow 
curves, in the limit of vanishing shear rates. Moreover, as 
discussed in more details in section IIII C[ the exponents 
found in Fig. [^b) and Fig. ^b), /3 = 4.0 and 5.75 for 
batch 1 and 2 respectively, only depend on the batch and 
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FIG. 3: (a) Shear rate 7 vs time for a = 32, 36, 37, 39, 41, 
42, 44, 45.5, 47, 50, 55, and 60 Pa (from bottom to top) for 
batch 2 and smooth BC. (b) Fluidization time 77 extracted 
from the second inflection point of y(t) (o) and from USV (•) 
vs the reduced shear stress a — o c with a c = 28.1 Pa. The 
gray line corresponds to the best fit 77 = B/{a — a c ) with 
B = (1.74 ± 0.04) ■ 10 10 and (3 = 5.75 ± 0.15. The inset shows 
the least square determination of a c (see text). 

not on the gap width or on the boundary conditions. We 
shall therefore conclude that the fluidization exponent is 
only function of the sample preparation, most probably 
through the final value of the pH of the microgel, as al- 
ready observed under imposed shear rate [ill ]. 



C. Influence of the geometry and carbopol 
concentration 

Creep experiments were performed in both rough and 
smooth plate-plate geometries for various mass concen- 



TABLE I: Evolution of the fluidization exponent /3 and pref- 
actor B with the carbopol concentration C. 



C (% wt.) 


B (Pa^.s) 
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0.5 


1.48-10 a 


3.4 ±0.1 
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17.44-10 5 
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FIG. 4: Fluidization time 77 extracted from the second in- 
flection point of "f(t) vs the reduced shear stress a — a c for 
various carbopol concentrations in a rough plate-plate geom- 
etry of gap 1 mm: [symbol, % wt., a c (Pa)]: (♦, 0.5, 17.3); 
(A, 0.7, 30.3); (•, 1, 30); (■, 3, 56.5). For the same batch 
of concentration 1 % wt., we test the influence of the bound- 
ary conditions and of the gap width: (symbol, gap, BC): (•, 
1 mm, rough); (o, 3 mm, rough); (®, 1 mm, smooth). Here, 
smooth BC correspond to glass plates. 



trations: C = 0.5,0.7,1 and 3 % wt. In this case, the 
samples are free of seeding glass spheres. The shear rate 
responses ^(t) (not shown) present the same characteris- 
tics as those of Fig. [^a) and Fig.[3fa). The corresponding 
fluidization times, extracted from the second inflection 
point of 7(t), are shown in Fig. [4] First, for the four con- 
centrations explored, we observe a power-law similar to 
the one found above in the Couette geometry. The expo- 
nent is an increasing function of the carbopol weight frac- 
tion C (see table |T|. Second, we checked on a sample of 
concentration C = 1 % wt. that the power law does not 
depend significantly on the gap width (1 < e < 3 mm) 
or on the BC. Moreover, since the roughness of the plate 
surfaces was changed either by using glass plates or by 
gluing sandpaper, we also infer that the chemical prop- 
erties of the surfaces do not play any noticeable role in 
the fluidization process. Thus, we conclude that for a 
given batch, the power-law behaviour of r/ vs a does not 
depend on the gap width or on the BC. 



IV. DISCUSSION 
A. Ruling out any fracture-like scenario 

In light of previous works, let us first justify our choice 
of a power law to fit the 77 vs a data. Indeed various 
other scalings have been proposed in the literature to 
account for stress-induced fluidization. Analogies with 
fracture in elastic solids [1] lead to 77 = r exp [(cro/cr) m ], 
with m = 1 for time-dependent rupture in glasses and 
m = 2 (to = 4 resp.) for nucleation of critical cracks 
that follow the Griffith's criterion in 2D (3D resp.) ge- 
ometries [29( . Activated processes have also been invoked 
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FIG. 5: Same data set as in Fig. [2Jb) plotted in various ways: (a) Fluidization time Tf vs the inverse of the shear stress a. The 
grey line corresponds to the best linear fit in a log-lin plot corresponding to Tf = to exp(<To/<7) with to = 11 ms and oo = 485 
Pa; (b) Fluidization time Tf vs the inverse of the shear stress squared a 2 . The grey line corresponds to the best linear fit in 
a log-lin plot corresponding to Tf = roexp [(cro/cr) 2 ] with to = 1 s and ao = 110 Pa; (c) Fluidization time Tf vs 1/cr 4 and 
(d) Tf vs a. The last two semilogarithmic plots allow us to rule out the corresponding behaviours, 77 = roexp [(cro/cr) 4 ] and 
Tf = tq exp(— a/ao) respectively. 



to support an exponential scaling ry = To exp(— a/ao) 
found in colloidal gels made of attractive carbon black 
particles [3(|. Fig. [5jc) andOJd) clearly show that the 
last two scalings cannot account for our carbopol data. 
Still fracture-like processes with m — 1 or m = 2 do fit 
our data correctly as found in [8] but the small range 
of a (less than one decade) does not allow for a definite 
conclusion derived from rheology alone. 

Here, however, one can go one step further since ul- 
trasonic echography provides crucial clues on the local 
yielding dynamics that rule out a fracture scenario. In- 
deed Fig. HI which is typical of all shear stresses in 
rough BC, shows that the backscattered pressure sig- 
nals do not present any discontinuity in the bulk from 
one incident pulse to another. This allows us to ex- 
clude the presence of fractures at least on scales above a 
few microns, in contrast with previous observations on a 
thixotropic organogel (see for instance Fig. 16 and Fig. 17 
in reference [3l|V Moreover, at short times (t < 60 s), the 
slopes of the echoes progressively decrease from the rotor 
to the stator, which is the signature of the homogeneous 
creep deformation seen in Fig. QJa). For t ~ 60 s, the 
whole sample suddenly stops as the material fails at the 
rotor (see the discontinuity at r ~ for t ~ 60 s) and 
subsequently undergoes total wall slip with a vanishing 
local shear rate. We interpret the complex shape of the 
echoes in the slip regime (see, e.g., the red box in Fig. [5]) 
as the consequence of elastic bulk deformations due to 
the large roughness (60 fim) of the moving wall. Indeed 



such erratic motions are not seen with smooth BC. Fi- 
nally the material also fails at the stator at t ~ 120 s 
(red arrow in Fig. [5]) so that it slips at both walls and 
undergoes a solid-body rotation, which is also typical of 
the slip regime in smooth BC. We can thus keep in mind 
that fracture planes similar to those found in solids can be 
ruled out and that carbopol microgels rather fluidize con- 
tinuously in a liquid-like fashion through the slow growth 
of a transient shear band. 



B. What can we learn from the creep deformation 
regime observed with rough BC? 

Let us now recall that at short times under applied 
shear stress and for rough boundary conditions [Fig. Q] 
and Fig. Ufa)], the shear rate follows a slow decay char- 
acterized by a power law j(t) oc i~ 2 / 3 known as the 
Andrade's law and often observed in solid materials 
[27l HH ■ Such a result may seem to contradict our previ- 
ous paragraph in which we have dismissed the idea of any 
fracture-like scenario or solid-like behaviour in the flu- 
idization process of carbopol microgels. Actually, there 
is no contradiction and we shall explain why in this sec- 
tion. 

In crystalline materials, creep deformation originates 
from the collective dynamics of dislocations (2(| which 
is intermittent at the mesoscopic scale (33j | but results 
in an average slow power-law relaxation of the mate- 
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rial at a macroscopic scale widely known as Andrade's 
creep [34]. The discrete dislocation dynamics model 
has revealed that sufficient -if not necessary- ingredients 
to observe Andrade's scaling law appear to be (i) long 
range aniso trop ic interactions and (ii) topological con- 
straints [lEISI!]- Dealing with (i), carbopol microgels are 
composed of highly crosslinked molecules and thus cer- 
tainly present long range interactions. Besides, the fact 
that we found a fluidization time that scales as a power 
law of the viscous stress could be a signature of such 
long range interactions, as suggested by numerical re- 
sults based on fiber bundle models [35j . Concerning (ii), 
carbopol microgels are amorphous materials and thus do 
not present any dislocations with well-defined topologi- 
cal constraints. Still, since carbopol microgels are made 
of soft jammed particles linked to each other by poly- 
mers, some topological constraints could rather apply to 
groups of particles presenting a low connectivity, so that 
both ingredients (i) and (ii) may be present in our soft 
system. 

A last relevant question remains the potential link be- 
tween the primary Andrade creep regime and the power- 
law scaling of the fluidization time: does this creep 
regime govern the behaviour of r/ or not? For rough 
BC, we observed such an initial creep regime and a flu- 
idization time decreasing as a power law of the viscous 
stress a — a c (Fig. [5]). For smooth BC on the other hand, 
there is no Andrade creep regime and the microgel enters 
immediately a total wall slip regime [Fig.|3[a)]. However 
we still observe a fluidization time that decreases as a 
power law of the viscous stress [Fig. |3Jb)]. Furthermore, 
the exponent does not seem to depend on the bound- 
ary conditions (Fig. [JJ. Finally, it can be checked in 
Fig. H] and Fig. [2fa) that the duration of the Andrade 
creep regime always remains at least one order of magni- 
tude smaller than the fluidization time. Altogether, these 
findings prove that for carbopol microgels, the Andrade 
creep regime certainly plays a negligible role on the flu- 
idization power law, contrary to what is seen for fiber 
composite materials where the primary creep regime is 
linked to the failure time of the samples [271] . This last 
point strongly distinguishes the fluidization of carbopol 
microgels from the failure process of solid materials. 



C. Comparison with experiments under controlled 
shear rate 

Under imposed shear rate, the fluidization scenario of a 
carbopol microgel is as follows [IJ HH : f° r rough bound- 
ary conditions and in a Couette cell geometry, the micro- 
gel first deforms elastically then undergoes plastic defor- 
mation and breaks at the rotor before entering a short 
regime of total wall slip. Then a shear band nucleates at 
the rotor and slowly invades the whole gap until a ho- 
mogeneous steady state is reached, with a linear velocity 
profile. The fluidization time of the carbopol, defined as 
the time to observe a linear velocity profile, scales as a 
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FIG. 6: Spatiotemporal diagram of the backscattered pressure 
signals coded in linear gray levels for a = 47 Pa applied at 
t = (batch 1, rough BC). Each horizontal line corresponds to 
ultrasonic echoes from the seeding glass spheres as a function 
of the spatial distance r to the rotor for an incident pulse sent 
at time t. Ultrasonic pulses are sent every 0.5 s. The slopes of 
the echoes in this diagram can be directly interpreted in terms 
of velocity, a vertical stripe corresponding to a vanishing local 
velocity. 



power law of the applied shear rate: Tf oc 7~ Q . ft is cru- 
cial to emphasize that the whole fluidization sequence is 
the same under controlled shear rate and under controlled 
shear stress and that in both cases the fluidization time 
decreases as a power law of the applied variable. The 
only difference between these two fluidization processes 
lies in the duration of the total wall slip regime, which 
is much shorter in the case of applied shear rate experi- 
ments. In the following, we make a connection between 
the two fluidization dynamics and the steady-state be- 
haviour. 



D. Link between transient and steady-state 
rheology 

1. From transient to steady-state rheology. 

The use of the power law = tj = B/(a — a c )P to 
describe stress-induced fluidization allows for a remark- 
able connection with independent experiments under im- 
posed shear rate and global steady-state rheology. In- 
deed, as recalled above, controlled shear rate fluidization 
on the same carbopol microgels was also shown to in- 
volve transient shear banding governed by a power-law 
behaviour: ri — Af-y , where a only depends on the 

batch preparation [TTJ] . Since rj and characterize 
the same physical process, we may crudely, but rather 
naturally, write that the two fluidization timescales are 
proportional: 

r^=XrP, (3) 



8 



160 
120 

cc 

S> 80 

D 

40 

\ Batch 1 ; rough BC (a) 

o' 

160 
120 

15" 

^ 80 
40 


0.01 0.1 1 10 100 

y(s" 1 ) 

FIG. 7: (a) Flow curve a vs 7 obtained by decreasing 7 from 100 to 0.01 s with a waiting time of 30 s per point for batch 1 
and rough BC. The gray line is the best fit by the HB model [Eq. fllj] with a c = 27.8 Pa, n = 0.53, and fj = 11.3 Pa-s~ n . The 
black dashed line is the HB model derived from fluidization times with a c = 28.3 Pa, n = 0.57, and fj — 9.4 Pa-s _n (see text). 
Inset: fluidization time Tj vs 7 in controlled shear rate experiments for batch 1 in rough BC, from Divoux et al. 11]. The 
line is = A/j a with A = 472 ± 11 and a = 2.30 ± 0.10. (b) Same as (a) for batch 2 and smooth BC. The best HB fit 
(gray line) performed on 0.2 < 7 < 100 s _1 yields a c = 30.4 Pa, n = 0.53, and fj = 10.3 Pa-s _n and fluidization times (black 
dashed line) lead to a c = 28.1 Pa, n = 0.51, and fj = 11.7 Pa-s~". The shape of this flow curve for 7 < 0.1 s _1 is typical of slip 
phenomena [36143S1 ] . Inset: for batch 2 in smooth BC with e = 0.5 mm, the best power-law fit of ri vs 7 yields A = 1375 ±46 
and a = 2.93 ± 0.15. 





where A is a dimensionless constant. This proportionality 
directly leads to: 

a = a c + fjj n , (4) 

with n = a If3 and fj = (B / XA) 1 ^ . In other words the 
HB behaviour that characterizes the steady- state rheol- 
ogy of our simple YSF is recovered from the two power- 
law behaviours of the fluidization times in independent 
transient experiments. Let us first test this result quan- 
titatively before trying to build it on more solid ground. 



2. Quantitative test. 

This link between local fluidization dynamics and 
global rheology is tested quantitatively in Fig. [7J for 
batches 1 and 2 and different BC. Going back to the ri 
vs 7 data (insets of Fig. [7]), we predict n — 2.3/4.0 = 0.57 
for batch 1 and rough BC and n = 2.93/5.75 = 0.51 for 
batch 2 and smooth BC. These exponents are in very 
good agreement with those extracted from the best HB 
fits of the flow curves, n = 0.53. It is also quite re- 
markable that the values of a c inferred from our fitting 



procedure of Tj (a c — 28.3 and 28.1 Pa, Fig. [2]) coincide 
with those of the best HB fits (a c = 27.8 and 30.4 Pa) 
to within 10 %. As a consequence the HB models in- 
ferred from the fluidization times are undistinguishable 
from the best HB fits of the steady-state rheology pro- 
vided A is left as a free parameter to achieve the correct 
values of fj from the prefactors A and B of the two power 
laws. We get A = 10.2 for batch 1 and 9.1 for batch 2, 
which means that fluidization under a given a is roughly 
ten times slower than fluidization under the correspond- 
ing 7. This difference most probably results from the 
fact that our YSF undergoes a long-lasting slip regime 
under stress while the slip regime is much shorter under 
imposed shear rate. 

3. Discussion. 

Let us now further justify writing that the two fluidiza- 
tion timescales and are proportional. First, we 
emphasize once again that both fluidization processes un- 
der imposed shear rate and under imposed shear stress 
are rigorously the same in that they present the exact 
same sequence of local events. Second, although the first 
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FIG. 8: Flow curves, shear stress a vs shear rate 7, obtained by decreasing 7 from 100 to 10~ 3 s _1 (o) and then increasing 
7 from 10 -3 to 100 s _1 (red line) for various waiting times per point: (a) t w — 2 s, (b) t w — 10 s, (c) t w = 30 s, and (d) 
t w = 70 s. Note that the hysteresis loop for 7 < 1 s _1 decreases for increasing waiting times. The total duration of the 
longest measurements shown in (d) is 2.1 ■ 10 4 s. Experiments performed in a plate-plate geometry (e = 1 mm) with rough BC 
(S = 46 (im) on a 1 % wt. carbopol microgel without any seeding glass spheres. 



two steps of the fluidization process, namely the homo- 
geneous strain and the total wall slip regimes, present 
a shorter duration in the case of applied shear rate, the 
fluidization times are of the same order of magnitude for 
both applied shear rate and applied shear stress. Third, 
a proportionality law is the simplest link one can think of 
between two timescales without artificially invoking an- 
other characteristic time whose origin and interpretation 
would be unclear. The fact that we can make any addi- 
tional timescale irrelevant will be confirmed in the next 
section. 

Now, our interpretation of the HB exponent as the ra- 
tio of two fluidization exponents allows us to revisit some 
recent observations. Indeed, quite similar results have 
been observed on acidic solutions of type I collagen by 
Gobeaux et at [39]. On the one hand, the steady-state 
rheology of type I collagen molecules is well described by 
a power law a — Kj n , with n ~ 0.27, i.e. a HB model 
with no yield stress. On the other hand, creep experi- 
ments have revealed that the fluidization occurs after a 
lag time t c that scales as t c oc l/er 1 /". This can also be 
read as t c cx l/(er — Ce) 1 /" with a c = 0. Using an origi- 
nal time stress superposition principle, the authors have 
shown that a and n are equal within error bars. Here, 
our approach allows us to predict that acidic solution of 
the same type I collagen should also fluidize under an ap- 
plied shear rate 7 after a lag time t c which should scale 
as t c oc 1/7^ with (3 = 1, so that if one imposes the lag 
time for both fluidization processes to be propotional, 
one would recover the steady-state rheology with n = a. 
Such a prediction remains to be experimentally verified 
and could reinforce our claim that the two fluidization 
timescales are proportional. 



Last but not least, in a recent work on bidimensional 
wet foams, Katgert et al. [40J also proposed an interpreta- 
tion of the HB exponent n. The authors made a connec- 
tion between the way the average drag force on a bubble 
scales with the velocity and the power-law behaviour of 
the viscous stress a — a c in the HB model. In this frame- 
work, n appears as a direct measure of the average forces 
at the bulk level taking into account both the local inter- 
bubble drag and the disorder induced by the flow. This 
is also compatible with our results, as we have observed 
that the fluidization exponents a and /3 are sensitive to 
the size distribution and/or the spatial organization of 
the soft particles that constitute the microgel, through 
the preparation protocol and the carbopol concentration 

(Fig. ED. 



E. Time-dependent effects and hysteresis cycles 

In this subsection, we provide the reader with further 
evidence showing that it would be artificial to invoke any 
other characteristic time or any hidden dynamical vari- 
able besides r/ to describe the fluidization process of car- 
bopol microgels. This point is strongly related to the fact 
that these microgels are simple YSF. 

Thixotropic and simple YSF are usually placed in two 
categories which preclude one another [2| . If satisfying at 
first sight, such a rough description remains qualitative 
and has resulted in defining useless sub-categories such 
as "unusual yield stress fluids" [4l[ . In order to overcome 
those difficulties, Coussot and Ovarlez [42[ recently pro- 
posed an interesting reunification of these two categories 
within a single theoretical framework by introducing the 
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ratio D of two timescales: a characteristic relaxation time 
Vo/Gq (built on the viscosity 770 and the elastic modulus 
Go of the fluid) and the restructuration time 9 of the sys- 
tem. A fluid with D = rjo/(Go6) close to 1 would be a 
simple YSF for which the restructuration time is indeed 
roughly equal to the relaxation time, whereas a fluid with 
D < 1 would correspond to a thixotropic material, con- 
comitantly presenting aging effects and restructuration 
over long durations. The key point of such a description 
is that one goes continuously from one type of fluid to 
the other simply by tuning the timescale ratio D. 

A simple way to probe the relevant timescales con- 
sists in performing successive decreasing and increasing 
ramps of controlled shear rates. Fixing the shear rate 
range (here, 7 m i n = 10~ 3 < 7 < 7 max = 10 2 s _1 ) and the 
number of experimental data points (here, 15 points per 
decade), the only control parameter is the waiting time 
per point t w spent at each imposed value of 7. In other 
words, i" 1 is the rate at which we scan the flow curve 
a(j). In Fig. [SI we report flow curves obtained with four 
different values of t w on a 1 % wt. carbopol microgel. 
As already briefly mentioned in [11] , for a given value of 
t w , we observe a slight hysteresis between decreasing and 
increasing shear-rate sweeps. Let us emphasize here that 
such an effect is also noticeable but not discussed in pre- 
vious literature [H, l2lj and that it is thus not particular 
to the type of carbopol that we are using. 

Repeating the shear rate sweep for different waiting 
times, we observe that the area A of the hysteresis loop, 
defined as 



A 



log 7m a, 



log 7n 



<r[log(7 ')]dpog(7 ')]. 



(5) 



decreases as power law of t w with an exponent 0.36 
[Fig. First note that this result is robust: the same 
trend is observed for carbopol microgels of three differ- 
ent mass concentrations (G = 1, 2, and 3 %) and, when 
rescaled by the elastic modulus Go of the microgel, A is 
roughly independent of G. Second, in light of the work by 
Coussot and Ovarlez [42| detailed above, this result is (i) 
compatible with a simple YSF, ruling out any thixotropic 
behaviour, and (ii) suggests that no other timescale than 
the fiuidization time 77 is necessary to describe the flow 
behaviour of carbopol microgels. Indeed, in the case of 
a simple YSF, for t w > 1 s, t w is generally large enough 
compared to the characteristic relaxation time of the fluid 
so that t w 3> 770/G0 ~ 0. Thus, the larger the imposed 
value of t w , the less the effects of the fluid relaxation 
and restructuration will be probed and so the smaller 
the area of the hysteresis loop. In other words, the larger 
t w , the more the carbopol microgel "forgets" about its 
shear history. 

In the case of a thixotropic YSF, the restructuration 
timescale would be of hundreds of seconds or more so that 
770 /Go <C 9 and t w would be of the same order of mag- 
nitude as 9. Thus, in the thixotropic case, at least two 
timescales, 9 and t w , would be involved in the material 
dynamics during shear rate sweeps and one would ex- 
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FIG. 9: Area A of the hysteresis between the decreasing and 
the increasing flow curve vs the waiting time per point t w for 
various carbopol weight fractions (symbol, % wt. carbopol): 
(•, 1 %); (J, 2 %); (□, 3 %). A decreases as a power law 
of the waiting time: A/Go = 0.41/tS; 36 , where Go is the 
elastic modulus of the microgel. Data obtained in a plate- 
plate geometry (e = 1 mm) with rough BC (8 = 46 /im) on 
carbopol microgels without any seeding glass spheres. 



pect a more complex behaviour of A vs t w than a simple 
decreasing function. In particular, it is anticipated that 
the hysteresis loop grows larger as long as t w < 9 and 
that it decreases (or at least saturates) for t w 3> 9 when 
the restructuration timescale is no longer relevant. An 
extensive study of the behaviour of hysteresis cycles for 
both simple and thixotropic YSF is currently underway 
to deeper test these ideas. In any case, to us, the results 
shown in Fig. [9] provide a strong confirmation that (i) 
carbopol microgels are simple YSF which exhibit negli- 
gible hysteresis when t w is large enough and (ii) no other 
timescale or hidden parameter is needed to describe the 
microgel rheology as the hysteresis can be accounted for 
only by a transient shear banding phenomenon which is 
fully described by the timescale 77. 



V. SUMMARY, OPEN QUESTIONS, AND 
OUTLOOK 

A. Summary 

We have performed a temporally and spatially resolved 
study of the stress-induced fiuidization of a simple yield 
stress fluid. The fiuidization is a four step process that 
successively involves Andrade-like creep deformation, a 
total wall slip regime, and a transient shear banding phe- 
nomenon that leads to a homogeneous flow in steady 
state. The time to reach a linear velocity profile is a 
robust decreasing power law of the applied shear stress 
which neither depends on the boundary conditions nor 
on the gap width, while the exponent is a function of the 
microgel microstructure. One of the key results of this 
article is that the exponent n in the HB model which de- 
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scribes the steady-state rheology naturally appears as the 
ratio a//3 of two fluidization exponents derived from in- 
dependent experiments under controlled stress and under 
controlled shear rate. To our knowledge, this provides for 
the first time a clear link between the transient regime of 
the fluidization process and the steady-state rheology. 

B. Open questions and outlook 

We would like to speculate that this last result is gen- 
eral for simple YSF and future experiments will focus 
on measuring and ri in emulsions and wet foams 
so as to extract the value of the exponent a and /3 and 
test their link with the steady-state rheology. Concerning 
Carbopol microgels, it would also be of valuable interest 
to unambiguously link the microscopic properties of the 
microgel, in particular the size of the microstructure, to 
the value of the fluidization exponents. 

The transient shear banding scenario, common to both 
applied shear stress and shear rate experiments, also re- 
mains to be characterized at a microscopic scale. In the 
case of traditional steady-state shear banding, the two 
flowing bands present two different microstructures [43| . 
In the case of wormlike micelle solutions for instance, 
the highly sheared band presents a nematic-like order 
whereas the micelles are more entangled in the weakly 
sheared band. Here, for the transient shear banding ob- 
served during the fluidization of carbopol microgels, one 
may wonder if there is any structural difference between 
the flowing band and the arrested region. 

Another puzzling issue comes up when one compares 
the fluidization laws of two different soft systems: car- 
bopol microgels whose fluidization time decreases as a 
power law of the viscous stress, and weakly attractive car- 
bon black gels whose fluidization time decreases exponen- 



tially with the applied stress [3(J ■ While the latter system 
is a fractal colloidal gel with a low volume fraction, car- 
bopol microgels are constituted of jammed swollen parti- 
cles. How and why does such a structural difference lead 
to different stress-induced fluidization law? Could one 
tune continuously the system properties to switch from 
one fluidization behaviour to the other? 

Finally, we wish to emphasize that it would be very 
interesting to compare the present experimental data on 
stress-induced fluidization to theoretical predictions. Un- 
fortunately, to the best of our knowledge, most recent 
theoretical works on shear banding in yield stress mate- 
rials have focused on stationary states only [42, 44] . One 
may think of using standard models for time-dependent 
materials based on structure-dependent kinetic equations 
[45M47| but, as discussed above, the restructuration ki- 
netics does not appear as a relevant ingredient for sim- 
ple YSF such as carbopol microgels. Two recent theo- 
retical papers 0, H^] based on the soft glassy rheology 
(SGR) model do address the transient regime but a de- 
tailed study of the fluidization times remains to be per- 
formed. Therefore, no analytical or numerical prediction 
is yet available for the fluidization behaviour observed in 
the present experimental work. Such a prediction would 
constitute a major step towards a full understanding of 
the yielding dynamics in simple YSF. 
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